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Abstract 

We report experimental measurements of the phase behavior of mixtures of thin (charged semi- 
flexible fd virus) and thick (/d-PEG created by covalently grafting poly- (ethylene glycol) to the 

o. 

surface of fd virus) rods. The /d-PEG are sterically stabilized but the fd virus are charged, thus 

pg ' varying the ionic strength of a mixture of fd and /d-PEG varies the effective diameter of the bare fd 

>' 
l/-) • rods, and the effective diameter ratio (d = Dfd-PEG/Dfd)- We examine the phase diagrams of the 

t> ; 

rod mixtures and find isotropic-nematic, isotropic-nematic-nematic and nematic-nematic coexist- 
ing phases with increasing concentration. In stark contrast to predictions from earlier theoretical 
work, we observe a nematic-nematic coexistence region bound by a lower critical point. Moreover, 
we show that a rescaled Onsager-type theory for binary hard rod mixtures qualitatively describes 
the observed phase behavior. 

PACS numbers: 64.70.Md, 61.30.St,61.30.Dk 
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The entropy driven phase transition of monodisperse suspensions of purely repulsive rods 
rom an isotropic to an aligned nematic phase has been extensively studied experimentally 
l|, |2|, theoretically J3, |4(, and computationally p. Binary mixtures of hard particles of 
different aspect ratios and/or shape can exhibit a much richer phase behavior. Quanti- 
tatively studying a binary hard-rod mixture is the first step towards understanding the 
impact of size polydispersity on the phase separation of concentrated suspensions of rodlike 
macromolecules. Length and diameter polydispersity are prevalent in suspensions of many 
synthetic and biological rodlike particles including F-actin, microtubules, and DNA. The- 
oretical studies of binary hard-rod mixtures predict that in addition to isotropic-nematic 
(I-N) coexistence, isotropic-nematic-nematic (I-N-N), isotropic-isotropic (I-I), and nematic- 
nematic (N-N) coexistence are possible when the length or diameter ratio of the particles 
is large enough p, LJ Lj lA llO, LLlI LL=1 ■ However, the topology of the theoretical phase di- 



agram of binary mixtures of rods, including the progression of the phase behavior from a 
totally miscible nematic to a demixed N-N state and the related existence of N-N critical 
points, is a subject of considerable debate as the predicted phase behavior of the more con- 
centrated phases (N-N and I-N-N) is extremely sensitive to the approximations employed 
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12J. lid- Il4| • The focus of past experimental studies has been on binary mixtures of rods 
of varying length where I-N, I-N-N and N-N coexistence have been measured |2lll5lllfil.ll7|. 
Because of polydispersity in the particle size, high solution viscosity and/or weak attrac- 
tions these measurements were constrained to low concentrations near the I-N transition 
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17J]. As a consequence, experimental studies of N-N phase behavior are rare. In this 
Letter, we present experimental measurements of the phase behavior of binary mixtures of 
rods of varying diameter ratio and equal length to very high (nematic) concentrations. These 
rods are strongly repulsive and highly monodisperse. By examining the phase diagrams well 
into the nematic region, we find a lower critical point in the N-N coexistence where previ- 
ously only an upper critical point had been predicted. We compare these unexpected results 
to the predictions of a scaled Onsager theory [12]. 

The theory describing the liquid crystal phase behavior of binary rod mixtures is closely 
related to the theory of orientational ordering developed by Onsager for monodisperse rods 
J2(]. Onsager showed that the decrease in orientational entropy accompanying the formation 
of an aligned nematic phase is more than compensated for by an increase in the free volume 
entropy of the rods. The I-N phase transition predicted for monodisperse rods is in excellent 



agreement with experiments |l[. By extending Onsager's theory for the free energy of hard 
rigid rods in the limit of the second virial coefficient (second virial theory, or SVT) to binary 



mixtures |3|, the phase behavior has 
diameter J9(, and different length 



3een determined for binary mixtures of rods of different 



J3, la |l4|. Onsager's second virial expansion of the free 
energy is, however, quantitatively accurate only for dilute suspensions of particles of large 
aspect ratio. For concentrated phase transitions (such as a N-N transition), both "end- 
effects" due to the finite size of the rods (typically neglected for Onsager-limit studies) and 
higher virial coefficients are significant, even in the limit of large aspect ratio. Therefore, 
to study the phase behavior of binary mixtures of rods of short, experimentally accessible, 
aspect ratios and to describe the phase behavior at high concentrations, we adopt a Parsons- 
Lee (P-L) free energy |l8L llfl| as applied to mixtures of rigid rods by Varga et al. |l2J |. The 
P-L approach approximates higher virial coefficients by interpolating between the Carnahan- 
Starling free energy for hard spheres and the Onsager free energy for long hard rods. The 
free energy is described by: 

^ = const + J2 Xl (\n(x lP ) + *(/,)) + |gl £ £ (%) 

8=1 HS t=l j=l 

Here p = (N% + N 2 )/V is the total number density, Xi and fi(fl) are the mole fraction and 
orientational distribution function of rods of type i, respectively, and a is the single particle 
orientational entropy |l2j |. In the second term, F^ is the residual free energy of a system 
of equivalent hard spheres of volume t>ns = (#1^1 + ^2^2) [12] , and vi is the volume of 
rod i. Furthermore, (vij) = J f XiXjVij(Q\, Q2)fi{^i)fj{^2)d^idQ2 is the full orientationally 
averaged excluded volume (vij) between rods i and j, with Q being the orientational unit 
vector (polar and azimuthal angles). The end-effects are included the excluded volume 
termQ. The theoretical phase diagrams presented here are calculated numerically with 
this free energy functional using techniques previously described |l2j]. 

Mixtures of two well characterized systems: the charged semiflexible fd virus Uj and fd 



virus irreversibly coated with the neutral polymer poly(ethylene glycol) (PEG) 2 (J, |2l( are 
investigated experimentally. The bare fd virus will serve as our thin rod, while the polymer 
coated fd ( fd-PEG) will be our thick rod. The fd virus is grown and purified as previously 



described 



22j, and it is characterized by its length L = 880 nm, diameter D = 6.6 nm, 



molecular weight M w = 1.64 x 10 7 g/mol, and persistence length, defined as the length over 
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FIG. 1: (Color Online) Phase separation of mixtures of fd and fd-PEG (20,000 M w PEG) at 110 
mM ionic strength {d = 3.7), pH 8.2. The phase diagram is presented in (a). The dark curves are 
guides to the eye representing the phase boundaries of the two nematic phases and the isotropic 
phase. The dashed lines indicate coexisting samples, and the three phase region is indicated by 
the grey triangle. The coexisting phases labeled in (a) as viewed under crossed polarizers are: 
(b) I-N coexistence; (c) highly fractionated I-N coexistence; (d) I-N-N three phase coexistence 
with a schematic representation of the partitioning of thick and thin rods shown in (e); (f) N-N 
demixing just above the triphasic region; and (g) highly concentrated N-N coexistence showing 
strong partitioning of the thick (yellow, /d-PEG) and thin (white, fd) rods. The isotropic phase 
appears dark and the nematic phase appears light due to its birefringence. The yellow color is 
due to the dye on the /d-PEG rods, (h) Differential interference contrast microscopy image of low 
density /d-PEG-rich nematic factoids in coexistence with high density fd-rich. nematic bulk phase. 

which the tangent vectors along a polymer are correlated, of P =2.2 jjlto. |l[. The fd rods 
also have a linear surface charge of 10 e~/nm at pH 8.2, thus to compare the phase behavior 
of fd with predictions for hard rods the electrostatic interactions are incorporated into an 
effective diameter D e g which is larger than the bare diameter and increases with decreasing 
ionic strength 2l|2j|. Over a wide range of ionic strengths, monodisperse suspensions of fd 
are known to undergo an isotropic to cholesteric (chiral nematic) phase transition which is 
accurately described by SVT for semiflexible hard rods of diameter D e g |l|. Because the 
free energy difference between the cholesteric and nematic phases is small, comparison with 



theories developed for nematic phases is appropriate 



24] 



The thick rods are formed by attaching an amino-reactive PEG, M w = 20000 g/mol 
(SSA-PEG20K, Shearwater Polymer Corp.), to the exposed amino termini of the viral coat 



jroteins. The dense polymer coating of approximately 200±30 PEG20K molecules per fd 



21 1 acts as a steric stabilizer above an ionic strength of 2 mM. In this regime the isotropic- 
nematic (cholesteric) phase transition of mono disperse suspensions of fd-PEG is independent 
of ionic strength J20|. The diameter -Dthick of the thick (virus+polymer) fd-PEG rods is 
thus calculated from the monodisperse isotropic coexistence concentration (q) using the 
theoretical relationship c ; [mg/ml] = 222/D thick [nm], as calculated from SVT adapted for 
semiflexible rods 20j, |25j. -Dthick varies from 25-40 nm depending on the /ci/SSA-PEG20K 
reaction conditions. Rods with different .Dthick were used separately. 

Samples are prepared at multiple virus compositions and concentrations to determine the 
phase behavior of mixtures of fd and fd-PEG. The virus suspensions are dialyzed against 20 
mM Tris-HCl buffer at pH 8.15 with NaCl added to change the ionic strength such that the 
diameter ratio D'thick/D'thin — d, where -Dthin = Des, varied from 1.1 < d < 3.7. The observed 
phase separation includes either an isotropic phase (I) coexisting with a nematic (N) phase, 
I-N-N three-phase coexistence, or N-N coexistence. The coexisting phases, as viewed under 



crossed polarizers, are depicted in Fig. [T] 
stable coexisting phases for such a system 



his confirms the theoretical predictions for the 



% |l2j . Phase separation between two isotropic 
phases (I-I) is not observed, presumably because d is not large enough. In Figs. ((It) and 
(0i) the /d-PEG-rich nematic phase can be seen to float above the fd-rich isotropic phase. 
Even though the volume fraction of rods is higher in the nematic phase, the mass density 
of the fd rich isotropic phase is greater than that of the fd-PEG rich nematic phase. The 
mass density difference arises in part because of the difference in single particle densities 
p = 1.35 g/ml for fd and p = 1.007 g/ml for fd-PEG. After equilibration, the concentrations 
of fd and fd-PEG in the coexisting phases are measured by absorption spectrophotometry. 
The optical density (A) of fd is A^ ^\ , for samples 1 cm thick. To independently measure 
the concentrations of thin and thick rods in the coexisting phases, fd-PEG is also labeled 
with the fluorescent dye fluorescein isothiocyanate (FITC) which has an optical density of 
^68000L/moi- Before bulk separation of coexisting phases occurs, phase separation on a micron 
length scale can easily be seen using fluorescence and/or polarization microscopy as shown 
in Fig. GJi- 

In Fig. El we present the experimental phase diagrams of the binary mixtures of fd and fd- 
PEG at three ionic strengths; corresponding to diameter ratio d ranging from 1.1 < d < 3.0. 
For large diameter ratios, d > 3.0, marked partitioning is seen in the I-N coexistence region, 
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FIG. 2: (Color Online) Phase diagrams of suspensions of fd-PEG (20,000 M w PEG) mixed with fd 
at (a) 60 mM, (b) 50 mM, (c) 5 niM ionic strength, pH 8.2. The legend for all three diagrams is in 
(c) . The dark curves are guides to the eye representing the phase boundaries of the two nematic 
phases and the isotropic phase. The dashed lines indicate coexisting samples, and open triangles 
are single phase nematic samples as determined by fluorescence microscopy. The three phase region 
is indicated by the grey triangle. 

as indicated by tie lines in the I-N coexistence region which do not pass through the origin. 
With increasing overall concentration, a triangular I-N-N three phase region and a N-N 
coexistence region are found (Fig. EK)- Almost complete partitioning of the thick (fd-PEG) 
and thin (fd) rods is observed in the N-N coexistence region. As d decreases, the I-N-N 
coexistence region becomes smaller, and below about d = 3.0 the triphasic region vanishes. 
Between d ~ 3.0 and d ~ 2.0 the N-N coexistence is still present at high concentrations, 
even in the absence of a well denned I-N-N triangle, suggesting a N-N region bounded by a 
lower critical point. Partitioning of the fd and fd-PEG between isotropic and nematic phases 
also decreases in the absence of an I-N-N region, as indicated by the shortening of the I-N 
tie lines which do not radiate from the origin from Fig. EJd to|2fc. For diameter ratios below 
d ~ 2.0 only a single nematic phase is observed. 

To shed light on the measured phase separation, we present the predictions for the phase 
behavior from P-L theory J12J in Fig. El Four distinct types of phase diagrams (indicated 
by the Greek symbols) are predicted for a mixture of rod-like particles of equal length as a 
function of the diameter ratio and the length of the thick rods. The qualitative evolution 
of the phase behavior as a function of d for L/D t hick — 24, which corresponds to the experi- 
mental aspect ratio of the thick rods, is displayed in Fig. Eb-d. The region of I-I coexistence 
(5) is not observed experimentally and is not discussed here. For small d, corresponding to 
region a, a transition from an isotropic to homogenous nematic phase is predicted at low 
concentrations and a N-N coexistence region bounded by a lower critical point is predicted 
at high concentrations (Fig. Eb)- Region a confirms that an I-N-N coexistence region is not 
required for the existence of a region of N-N coexistence, as experimentally observed (Fig. 
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12b)- As d increases within region a the N-N lower critical point moves to lower rod concen- 
trations. Region (3 is characterized by a phase diagram which exhibits an I- N-N coexistence 
region capped by a region of N-N coexistence bounded by an upper critical point, with an 
additional region of N-N coexistence bounded by a lower critical point (Fig. |3fc). Upon 
increasing d, these two regions of N-N coexistence coalesce to form a single N-N coexistence 
region (region 7, Fig. EJi). 

The phase diagram for region 7 is qualitatively similar to that predicted by SVT [9| 
for d > 4.3, and is experimentally observed for d > 3 (Fig. |2K) • The phase diagrams 
characterizing regions a and ft however, are q nalitatively different from those pred.eted in 
the SVT (L/D — * 00) |9| because they predict a region of N-N coexistence bounded by a 
lower critical point. Only an N-N coexistence region bounded by an upper critical point is 
predicted in the SVT for 3.8 < d < 4.3 J3|. However, it is precisely this N-N coexistence 
region bounded by a lower critical point which is measured in the fd+fd-PEG solutions (Fig. 
I2t>) for 3 > d > 2. Simply by extrapolating between the Onsager limit (L/D — > 00) and 



the Carnahan-Starling hard sphere limit (L/D — > 1) 18], the P-L scaling |12J qualitatively 
reproduces the experimental phase diagram for 3 > d > 2 while SVT alone does not. 

Comparing the evolution of experimental phase behavior with d for our long semiflexible 
rods with the P-L predictions we find that it qualitatively follows the phase behavior pre- 
dicted for short rigid rods, L/D t hick < 7. In this case, region (3 is bypassed and the N-N 
region bounded by the lower critical point coalesces directly with the I-N region creating an 
I-N-N coexistence region. It has been shown in simulations that the excluded volume of a 



3. 



Thus 



flexible rod is equivalent to the excluded volume of a shorter but thicker rigid rod 
we expect long flexible rods to exhibit a phase behavior similar to that predicted for shorter 
rigid rods, as observed. Additionally, we observe that the experimental I-N-N coexistence 
is stable to much lower diameter ratios, d ~ 3, than predicted. We suggest that this is 
because the thin-thick rod interactions are non-additive. Because -D t hi n is defined by elec- 
trostatic repulsion, nothing prevents the polymer coated virus (thick rods) from touching 
the bare surface of the fd (thin rods). Subsequently, the thin- rod interaction diameter is ef- 
fectively smaller (d is effectively larger). Theoretically, one of the challenges that remains is 
to incorporate non-additivity and flexibility into theories for the binary rod phase behavior. 
In conclusion, we have studied the phase behavior of the first experimental system of 
thick and thin rods. Onsager's SVT qualitatively reproduces the main features seen in our 
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FIG. 3: (a) Characterization of the binary phase behavior of rods of equal length obtained with 
P-L scaling of the Onsager free energy as a function of the inverse diameter ratio 1/d and L/Di^^. 
I-N-N coexistence is found below the uppermost line (regions /3 and 7). Coalescence of the upper 
and lower N-N critical points occurs at the dashed line (region 7), and I-I-N coexistence becomes 
stable below the lower line (region 5). For L/£>thick < 7> the lower critical point merges with 
the I-N coexistence region at the uppermost solid line. The phase diagrams are calculated for 
L/Dthick = 24 with d = 4.4 (b), d = 5 (c), and d = 5.2 (d). Phase diagrams are presented as a 
function of reduced concentration c, = ViNi/V . Three phase coexistence is indicated by the grey 
triangle. 



experimental binary rod mixtures at large diameter ratios, but does not accurately capture 
the evolution of the phase behavior from a totally miscible nematic to a demixed nematic- 
nematic state. By incorporating the higher virial coefficients in an approximate fashion with 
the Parsons-Lee scaling of the free energy for hard rods |l2|], we proposed an evolution of the 
binary rod phase diagram with diameter ratio which qualitatively describes the experimental 
findings. Our experimental and theoretical results show that an I-N-N coexistence region is 
not required for the existence of a region of N-N coexistence in contrast to past predictions. 
However, the N-N upper critical point, which is predicted for very long rods in both the SVT 
and Parsons-Lee theory has not yet been observed experimentally; further experimental or 
computational studies of binary mixtures of longer, more rigid rods, may reveal this upper 
critical point. 
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